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Abstract 

We have searched for lepton family number violating decays of the r lepton 
using final states with an electron or a muon and one or more ir° or rj mesons 
but no neutrinos. The data used in the search were collected with the CLEO 
II detector at the Cornell Electron Storage Ring (CESR) and correspond to 
an integrated luminosity of 4.68 fb~ l . No evidence for signals was found, 
resulting in much improved limits on the branching fractions for the one- 
meson modes and the first upper limits for the two-meson modes. 
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In physics all fundamental conservation laws have associated symmetries. Lepton flavor 
conservation is an experimentally observed phenomenon with no associated symmetry in the 
Standard Model. Lepton flavor violation is expected in many extensions of the Standard 
Model such as lepto-quark, supersymmetry, superstring, and left-right symmetric models, 
and models that include heavy neutral leptons. Gonzalez- Garcia and Valle have calcu- 
lated 1], in a model with Dirac heavy neutral leptons, the branching fractions for r decay 
into one lepton plus a photon or a 7r° or 77 meson. The branching fractions depend on the 
heavy neutral lepton masses and mixings. Given the constraints from other measurements, 
the branching fraction for r~ — > £~tt° may still be as large as 10 -6 for neutral lepton 
masses above a few TeV/c 2 and is higher than that for the radiative decay t~ — > £~ / y. 
Using a Grand Unified Theory (GUT) and superstring inspired model with heavy neutral 
leptons, Ilakovoc and collaborators 0] have calculated the branching fractions for r decay 
into one lepton plus one or two mesons. The branching fractions depend on the masses of 
the Majorana neutrinos and the mixings between the heavy and light neutrinos. The decay 
t~ — > £~7t° may have a branching fraction as large as 10~ 6 . The previous upper limits j|] 
on the branching fractions for the decays into one lepton and a tt° or 77 meson are of the 
order of 10 -4 — 10~ 5 and for the decays into one lepton and two charged ir mesons are of the 
order of 10~ 6 . There are no published results for the decays into one lepton and two neutral 
mesons (7r°7r°, rjr] or 7r°77). The CLEO II experiment with its large sample of r events may 
have the sensitivity to observe the lepton flavor violating decays. In this Letter, we present 
the result of a search for the decays into one lepton and one or two 7r° or 77 mesons. 

The data used in this analysis were collected with the CLEO II detector from e + e~ 
collisions at the Cornell Electron Storage Ring (CESR) at a center-of-mass energy ~ 10.6 
GeV. The total integrated luminosity of the data sample is 4.68 fb _1 , corresponding to the 
production of N TT = 4.26 x 10 6 t + t~ events. CLEO II is a general purpose spectrometer || 
with excellent charged particle and shower energy detection. The momenta and specific 
ionization (dE/dx) of charged particles are measured with three cylindrical drift chambers 
between 5 and 90 cm from the e + e~ interaction point, with a total of 67 layers. These 
are surrounded by a scintillation time-of-flight system and a CsI(Tl) calorimeter with 7800 
crystals. These detector systems are installed inside a superconducting solenoidal magnet 
(1.5 T), surrounded by an iron return yoke instrumented with proportional tube chambers 
for muon identification. 

The r + r~ candidate events must contain exactly two oppositely charged tracks. To 
suppress beam-gas events, the distance of closest approach of each track to the interaction 
point must be within 0.5 cm transverse to the beam and 5 cm along the beam direction. 
We divide each event into two hemispheres (signal and tag), each containing one charged 
track, by the plane perpendicular to the thrust axis, which is calculated using both charged 
tracks and photons. The total invariant mass of the tag hemisphere must be less than the r 
mass (M T = 1.777 GeV/c 2 ) [|]]. Because there is no neutrino in the signal hemisphere while 
there is at least one neutrino undetected in the tag hemisphere, the missing momentum of 
the event must point toward the tag hemisphere. To suppress the background from radiative 
Bhabhas and fi pairs, the direction of the missing momentum is further required to satisfy 
\cos9 m i SS i n g\ < 0.90, where 9 is the polar angle with respect to the beam. To reject the 
background from two-photon interactions, we require the magnitude of the net transverse 
momentum vector of each event to be greater than 300 MeV/c. 
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The signal hemisphere must contain an electron or a muon. The electron candidate must 
have a shower energy to momentum ratio in the range, 0.8 < E/p < 1.1, and have specific 
ionization loss within three standard deviations of that expected for an electron. The muon 
candidate must penetrate more than three absorption lengths of material. 

Photon candidates are defined as energy clusters in the calorimeter of at least 60 MeV in 
the barrel region (|cos#| < 0.80) or 100 MeV in the endcap region (0.80 < \cos9\ < 0.95). We 
reconstruct 7r° and r] mesons using the 77 decay channel. For the decays involving one meson, 
both photons must be in the barrel. For the decays involving two mesons, at least one photon 
from each meson must be in the barrel. We further require every photon to be separated from 
the projection of any charged track by at least 30 cm unless its energy is greater than 300 
MeV. There is no explicit cut on the maximum number of photons in the signal hemisphere in 
order to maintain a high detection efficiency while minimizing the dependence on the Monte 
Carlo simulation of electromagnetic showers. The signal hemisphere may contain photons 
not used in the 7r°/?7 reconstruction. However, photon candidates with energy greater than 
300 MeV or with a lateral shower profile consistent with that expected for a real photon 
must be used in the reconstruction. For the one-meson (two-meson) mode, events with more 
than two (four) such photons in the signal hemisphere are rejected. 

To search for neutrinoless r decays, we select r candidates with invariant mass M and 
total energy E within the ranges, 

— 250 < AE = E — < 150 MeV, 

-80 < AM = M - M T < 60 MeV/c 2 . 

These requirements correspond approximately to three standard deviation limits, according 
to the Monte Carlo simulations (see below). We then look for 7r° and 77 candidates using the 
77 invariant mass spectrum. The mass spectrum is expressed in standard deviations from 
the nominal 7r° or rj mass 0], 

where <r 77 is the mass resolution calculated from the energy and angular resolution of each 
photon. The S" 77 distributions may have multiple entries due to different combinations of 
photons in the 7r°/r) reconstruction. The signal region is defined as —3 < S* 77 < 2 while the 
sideband regions are defined as —10 < S' 77 < — 5 and 4 < 5 77 < 9. Two signal events satisfy 
these selection criteria (see Table |). The AE vs. AM and 5" 77 distributions of these events 
are shown in Fig. |]. 

The detection efficiencies (e) are estimated using a Monte Carlo simulation. In the Monte 
Carlo, one r lepton decays according to a two-body (three-body) phase space distribution 
for the one-meson (two-meson) mode of interest and the other r lepton decays generically 
according to the KORALB r event generator || . The detector response is simulated using the 
GEANT program [^] . The background from generic r decays is estimated using the KORALB 
Monte Carlo and the background from hadronic events is estimated using the Lund Monte 
Carlo ||. The detection efficiencies and background estimates are summarized in Table Q. 
Also listed are background estimates based on the numbers of events in the sideband (N s b) 
and corner-band (N c b), if appropriate, regions in the 5 77 distribution, Nb g = ^N s b — \N c b ) 
assuming a linear background distribution. In the data, N c b is measured to be zero for all 
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TABLE I. Summary of detection efficiencies, signal, backgrounds, 90% C.L. upper limits on 
the signal (see text) and branching fractions. 



Mode 


e(%) 


N ob 


e-TT° 


8.71 




n 


8.08 






9.97 






8.49 




e-7r 7T 


5.12 




/i-7r vr 


3.71 


1 


e~rjrj 


6.09 






3.48 




e~ir°i] 


5.52 


1 




3.73 





Ao 


Ag 


Ap 


A 


6(10 b ) 


2.30 


2.32 


2.30 


2.32 


3.72 


2.30 


2.32 


2.30 


2.32 


4.01 


2.30 


2.32 


2.30 


2.32 


8.19 


2.30 


2.32 


2.30 


2.32 


9.62 


2.30 


2.34 


2.30 


2.34 


6.47 


3.89 


3.97 


3.61 


3.68 


14.0 


2.30 


2.34 


2.30 


2.34 
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2.30 


2.34 


2.30 


2.34 


60.2 


3.89 


3.96 


3.61 


3.68 


23.8 


2.30 


2.34 


2.30 


2.34 


22.4 



0.5 



0.5 



0.4 



0.2 0.4 



six two-meson decays. Because of the paucity of events, rather than comparing the number 
of events observed to the expected background in each individual mode, we will sum over 
all the modes for comparison. The two events observed is somewhat higher than the 0.5 
events expected from the sideband technique but consistent with the 1.5 events estimated by 
the Monte Carlo simulations. There is therefore no evidence for a signal. The background 
estimated using the sideband technique is used to compute the upper limit on the signal. 
The upper limit on the branching fraction is related to the upper limit A on the signal 

by 



where B\ is the inclusive 1-prong branching fraction B n o (B v ) is the branching fraction || 
for 7T° — > 77 (77 — > 77), and m (n) is the number of 7r° (77) mesons in the final state. The 90% 
confidence level upper limits on the signal are summarized in Table |. We calculate the upper 
limit A using a Monte Carlo technique, which incorporates both the Poisson statistics of the 
signal and the systematic errors. The systematic errors include the statistical uncertainty 
in the background estimate due to limited statistics in the sideband (and corner-band if 
appropriate) regions. This statistical uncertainty is incorporated using Poisson statistics. All 
other sources of systematic errors are incorporated using Gaussian statistics. These include 
the uncertainties in the t + t~ cross section (1.0%), luminosity (1.0%), track reconstruction 
(3.0%), lepton identification (1.5% for e and 4.0% for 71), 7r° or 77 meson reconstruction (5.0% 
per meson), branching fraction of 77 — > 77 (0.8%) ||, and detection efficiencies due to limited 
Monte Carlo statistics (2-3% for the one-meson modes and 3-4% for the two-meson modes). 
These uncertainties are added in quadrature in computing A. For comparison, we also list 
the upper limits A , Xq an d Ap. A is calculated using only Poisson statistics for the signal, 
Ag includes all the systematic errors except the statistical uncertainty in the background 
estimate, and Ap includes only the statistical error in the background estimate. Ag is larger 



6 



than Ao as expected. However, Xp is smaller than Ao when the observed number of events 
is non-zero and the estimated background is zero. This is not unexpected because in the 
calculation of Xp we allow for the estimated zero background to fluctuate up, in contrast to 
Ao in which the background is estimated to be zero with no uncertainty. 

The upper limits on the branching fractions for the modes involving one neutral meson 
are significantly more stringent than the published results [|J. There are no previous limits 
for the modes involving two neutral mesons. The limits for the 7r°7r° modes are comparable 
with the limits for the 7r + 7r~ modes Q. In the model of Gonzalez- Garcia and Valle 0, the 
limit on t~ — ► e~7r° extends the heavy lepton mass vs. mixing region previously excluded 
from other measurements: neutral leptons with mass greater than 6-10 TeV/c 2 for mixing 
with the third generation in the range 0.05-0.02 are now excluded ||. 

We gratefully acknowledge the effort of the CESR staff in providing us with excellent 
luminosity and running conditions. This work was supported by the National Science Foun- 
dation, the U.S. Department of Energy, the Heisenberg Foundation, the Alexander von 
Humboldt Stiftung, Research Corporation, the Natural Sciences and Engineering Research 
Council of Canada, and the A. P. Sloan Foundation. 
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FIG. 1. AE vs. AM and S" 77 vs. S 77 distributions of t — > fi 7r^7r" (a-b) and t — > e tt^t] (c-d) 
candidates in the data and signal Monte Carlo (open circles) samples. The 5 77 vs. S 77 distribution 
is for the center box in the AE vs. AM plane. The signal, sideband, and corner-band regions in 
the S-p, vs. S 77 plane are indicated by the 9 boxes. The size of the circles is proportional to the 
number of entries. The scale for the signal Monte Carlo event distributions is arbitrary. 
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